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Abstract. We describe a sample of thermally emitting neutron stars discovered in the ROSAT All-Sky Survey. We discuss the 
basic observational properties of these objects and conclude that they are nearby, middle-aged pulsars with moderate magnetic 
fields that we see through their cooling radiation. While these objects are potentially very useful as probes of matter at very 
high densities and magnetic fields, our lack of understanding of their surface emission limits their current utility. We discuss 
this and other outstanding problems: the spectral evolution of one sources and the relation of this population to the overall 
pulsar population. 
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INTRODUCTION 

While almost 2000 isolated neutron stars have now been 
discovered as radio pulsars, the total number in the 
Galaxy is much larger. Radio pulsars emit pulsations for 
~ 10 yr and are visible due to radio beams that subtend 
1-10% of the sky, so the total number of neutron stars 
of all ages just in the local region of the Galaxy (where 
radio pulsars are detectable) should be > 10 6 [e.g.,|H|2l. 

Are those neutron stars not detectable as radio pulsars 
objects invisible, or is there some chance of observing 
them (there is always the exception of Geminga, where 
we see no radio emission but is otherwise a standard 
pulsar)? For years astronomers have proposed that a large 
fraction of these objects would be visible through one 
of two mechanisms: accretion JH-ID] or cooling Jg]. The 
first mechanism could revive old, dead pulsars, while the 
second would primarily work for younger sources. Both 
mechanisms, however, make the neutron stars visible in 
the soft X-ray regime, not in the radio regime that had 
dominated the study of neutron stars. 

These neutron stars should be identifiable by J3]: 

1 . Largely thermal emission peaking in the soft X-ray 
or far-UV band, requiring small hydrogen column 
densities to remain visible 
The absence of bright optical counterparts 
Significant (> 0.1 arcsec yr -1 ) proper motions 
Preferred locations in the Galactic plane 



The first two criteria relate to the spectra of the neu- 
tron stars, and serve to rule out the active galaxies and 
stars that dominate X-ray surveys [7i\. The third crite- 
rion reflects the proximity of the sources (with maxi- 
mum distances of ~ 1 kpc) and the large space velocities 
of known neutron stars [2, 8]. The final criterion comes 



from the Galactic nature of the sources, and is similar to 
the distribution of radio pulsars. 



The Legacy of ROSAT 

While some had anticipated up to 5,000 objects dis- 
covered through soft X-ray surveys J3, 0], the ROSAT 
All-Sky Survey (RASS) discovered just over half a 
dozerQ (as of 2007) nearby cooling neutron stars (see 
Tab. [TJ. These neutron stars are then all the more valu- 
able because of their rarity. 

The first such source to be discovered was 
RX J1856.5-3754 (hereafter RX J1856). It was origi- 
nally identified serendipitously as a soft, bright X-ray 
source with no obvious optical counterpart JT2I1 . Its 
location in front of the R CrA molecular cloud meant 
that it had to be nearby (< 200 pc IU3I1 ) and hence small 
(not a white dwarf or anything larger) — otherwise the 
X-ray emission would have been absorbed^. Confirma- 
tion of its nature came with the discovery of a very faint 
(niB « 25.8 mag), blue optical counterpart ifnll . 

Since then, six other similar sources have been iden- 
tified with considerable effort 11371 142M46I1 . While a va- 
riety of names and acronyms exist for these objects, we 
call them simply "Isolated Neutron Stars" (INS). Iden- 
tification of additional sources that may still be present 
l47ll in the ROSAT Bright Sources Catalog (contain- 
ing w 18000 sources with > 0.05 counts s _1 with the 



1 The difference is largely from poor assumptions about the pulsar 
velocity distribution |6, 9] and the effects of magnetic fields [ HI [Till . 

2 As it turned out, this argument was false, as observations of more 
distant stars did not have high extinction flUl . Nonetheless, RX J 1856 
is quite close fl5l[l6ll . 



TABLE 1. Observed Properties of the Seven Isolated Neutron Stars 
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* We give the spin-period, period derivative, and rms pulsed fraction. 

* We give the hydrogen column density, blackbody temperature, XMM-Newton EPIC-pn count-rate, central energy of any absorption features in 
the X-ray spectrum, and B-band magnitude. We also give energies of any secondary lines, if known. Note that all spectral estimates are covariant, 
and that the details of the X-ray absorption depend on the specific model. 

** We give the proper motion and distance, using the parallax if known, else that inferred from Posselt et al. [41]. 

* The spectral parameters are average quantities. 
§ Inferred from STIS 50CCD. 

' Based on Motch et al., these proceedings. 
" Inferred from V and r 1 bands. 



Position-Sensitive Proportional Counter, or PSPC; JH]) 
is extremely difficult given the poor positional accuracy 
of the PSPC. We will return to this later. 



OBSERVED PROPERTIES 

We summarize the properties of the 7 confirmed INS 
in Table [TJ also see previous reviews such as Haberl 
l29ll . that give additional observational details. The basic 
properties largely bear out the expectations of Treves 
and Colpi [4] but with much smaller numbers: they are 
soft X-ray sources with faint optical counterparts and 
significant proper motions. However, there are too few 
of them to see a concentration in the Galactic plane, and 
they may actually reflect a more local population that was 
born in the Gould Belt JH]. 

The X-ray spectra of the INS are reasonably close to 
soft black bodies with kT in the range of 40 to 100 eV 
lfl8tl attenuated by a small amount of interstellar absorp- 
tion (Fig.[TJ). There is no sign of any non-thermal X-ray 
power-law such as those seen in the spectra of most radio 
pulsars. Most of the INS, however, have spectra that ap- 
pear with si gnifi cant absorption features at low energies 
(Fig. [D Hall USE! S3- Typically, the absorption is 
modeled as one or more broad absorption lines, although 
the line shapes are complex and phase-dependent, and as 
we obtain increasing amounts of data our initial fits are 
no longer sufficient. All but one of the INS have hints 
of absorption in their spectra, and several 11291 13611 may 
even have absorption at multiple energies. Posselt et al. 
l4lll used the measured column densities, along with a 
map of the local interstellar medium, to infer distances to 



the INS. They found the objects were nearby (< 500 pc, 
generally), although we note that the it can be difficult 
to establish a reliable column density when the intrinsic 
form of the spectrum is not well known [e.g.. l50ll . 

All but one of the INS show gentle, largely sinusoidal 
X-ray pulsations, although some may be double-peaked 
i35il . The periods are all tightly clustered (compared to 
radio pulsars) between 3 and lis (Fig.[2j, and the pulsed 
fractions vary from ~ 1% to almost 20%. Through re- 
peated X-ray observations we have been able to estab- 
lish reliable, phase-connected timing solutions for two 
of the INS, RX J0720.4-3 125 (hereafter RX J0720) and 
RX 1308.6+2127 (hereafter RX 1308), and we find spin- 
downs of ~ 10~ 13 s s , With the period derivatives, we 
can calculate the usual pulsar quantities l5lil : dipolar 
magnetic field Bdip = 2.4 and 3.4 x 10 13 G, character- 
istic age T —1.9 and 1 .5 Myr, and spin-down energy loss 
rate E = 4.7 and 4.0 x 10 30 erg s" 1 for RX J0720 and 
RX 1308, respectively HQ. 

In the optical and ultraviolet, deep searches have re- 
vealed confirmed counterparts to four of the INS lfl7[l27i 
HHHHIHl], an d a possible counterpart to a fifth 1 4011 . In 
keeping with the identification of these objects as neu- 
tron stars, the optical counterparts are quite faint (far 
fainter than any possible companion might be), with X- 
ray to optical flux ratios of ~ 10 4 . However, the opti- 
cal fluxes generally lie a factor of ~ 10 above the ex- 
trapolation of the X-ray spectrum (Fig.[TJ): the so called 
"optical excess." In the best studied case, RX J 1856, 
the optical/UV spectrum is consistent with the slope of 
a Rayleigh-Jeans tail: F v °= v 2 014II . In the other cases 
we do not have nearly as much data and the inferences 
are consequently less certain. For RX J0720, the opti- 




FIGURE 1. Spectral energy distributions for RX J 1856 (left), RX J0720 (center; with the X-ray spectrum from before the 
appearance of an absorption feature), and RX J 1605 (right) . For the first two, the X-ray points are from Chandra/LETG spectra 
fl8Ll2lll . and for RX J 1605 the X-ray points are from XMM-Newton/RGS [28] The optical data are from a combination of ground- 
based and space-based observations [ 14, 27], The solid curves represent the best-fit black-body models to the X-ray data (with an 
absorption feature in the case of RX J1605); the dotted curves are the same model without interstellar extinction. 



cal/UV spectrum is close to the slope of a Rayleigh-Jeans 
tail, but there are indications that it deviates lull l22tl . 
For RX J1605.3+3249 (hereafter RX J1605), the data 
are even sparser, but it might have even larger deviations 
from a Rayleigh-Jeans tail l30ll . 

An elongated Ha nebula surrounding RX J 1856 was 
found by van Kerkwijk and Kulkarni 15411 . The nature of 
this nebula is not entirely clear, but it is likely a bowshock 
formed by the interaction between the neutron star's en- 
ergetic wind and the interstellar medium through which 
it is traveling supersonic ally [e.g.,|55t]. In this model, we 
can infer E for RX J 1856 even though P has not been 
measured, and find E > 1.2 x 10 33 c/| 60 erg s _1 (where 
the distance is 16(Mi60 pc). Motch et al. lf30h discov- 
ered some evidence of Ha emission from the position 
of RX J1605, but it has yet to be confirmed. 

Proper motions have now been measured for four of 
the INS iH E2 El and Motch et al., these pro- 
ceedings]. The implied transverse velocities are gener- 
ally ~ 200 km s~ 1 , consistent with the pulsar population. 
Tracing the proper motions back, one finds that these INS 
are consistent with being born in nearby OB associations 
~ 1 Myr ago 15611 . We have also been able to measure 
parallaxes to two of the INS. The distances are largely 
consistent with those estimated by Posselt et al. Ir4~lll . 

In the radio, there are no confirmed detections of any 
INS as a steady radio source or via pulsations 1I2II 1271 
l32i l57ll . Motivated by the similar location in the P — P 
diagram of the INS and the Rotating RAdio Transients 
(RRATs; 15811 ). there have also been searches for spo- 
radic ratio bursts from the INS, and these too have been 
negative [see Burgay, Kondratiev in these proceedings]. 
Similarly, searches for counterparts in the near-infrared 
(where one might see emission from a faint companion 
or an accretion disk) have not been successful 05911 . 



INFERENCES & ENERGY SOURCES 

From the general properties of the INS discussed above, 
it seems clear that the INS are reasonably young (~ 
1 Myr) isolated neutron stars. The X-ray emission that 
we see is plausibly thermal: the values of E that we 
measure or infer are all quite low, compared to X-ray 
luminosities of - 3 x 10 32 c/| 60 erg s" 1 (for RX J0720 
12111 . where we normalize to a distance of 360o?36o pc). 
In contrast to radio pulsars (and other rotation-powered 
objects such as Geminga), which have Lx ~ 10^E l60l 
16111 and where rotational energy powers the majority of 
the X-ray emission (some may be thermal), RX J0720 
has Lx/E ~ 60. Therefore, for RX J0720 rotation cannot 
power the X-ray emission, and we must resort to other 
mechanisms. While the value of E inferred for RX J 1856 
from the Ha nebula is of a different order than those for 
RX J0720 and RX 1308, it is still much less than what 
would be necessary to contribute significantly to the X- 
ray emission. However, unlike for RX J0720 where the 
age inferred from timing is close to the kinematic age, 
using the inferred E for RX J 1856, its spin period, and 
kinematic age it is hard to form a consistent picture, 
something we are investigating. 

The remaining energy sources that we must consider 
are accretion, magnetic fields, and thermal energy, all of 
which are seen in other types of neutron stars. 

Models involving accretion from the ISM for the INS 
[e.g., IH] have essentially been ruled out as energeti- 
cally important: the high velocities (~ 200 km s -1 ) in- 
ferred for RX J 1856 and RX J0720 from their proper 
motions and parallaxes Il5l l24ll make Bondi-Hoyle ac- 
cretion (M °= v~ 3 ) very improbable, especially with the 
ISM density for RX J 1856 inferred from the Ha neb- 
ula H54fl . We can then examine accretion from a fossil 
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FIGURE 2. A portion of the P-P diagram, showing radio 
pulsars (points; 1 65]), magnetars (diamonds), and the newly- 
discovered rotating radio transients (RRATs 15811 ; stars and 
arrows at bottom); selected objects are labeled, and the pulsars 
with > 0.05 s in the RASS are shown as the squares. Also 
shown are the INS: RX J0720 and RX J 1308 by filled circles, 
and the others by arrows at the top (since P is unknown). 
The diagonal lines show loci of constant dipole magnetic field 
(solid) and spin-down age (dotted), while the dot-dashed line is 
an approximate "death-line" for pulsar activity. The similarity 
in spin properties between the INS and some of the RRATs was 
noted by McLaughlin et al. 1 58] and Popov et al. 1 66]. 

fall-back disk JUB- Timing observations of RX J0720 
limit the torques from a disk such that < 1% of the ob- 
served X-ray luminosity can come from frictional heat- 
ing, and having accretion penetrate the centrifugal bar- 
rier to land on the surface would require an unreasonably 
large accretion rate of M ~ 2 x 10 17 g s _1 [also see l59tl . 

Excluding accretion as an energy source, we are left 
with residual thermal energy and magnetic fields. All 
neutron stars should have thermal energy powering emis- 
sion at some level, but only those with magnetic fields 
high enough to decay (> 10 14 G, generally) will have 
significant magnetic luminosity ll67ll . 

Heyl and Kulkarni l67ll proposed a model where 
RX J0720 was heated by magnetic field decay. This 
would help explain the similarity between the spin- 
period of RX J0720 and those of magnetars l68Tl . and 
it would also explain why objects like RX J0720 could 
be overrepresented in a local sample. However, timing 
observations of RX J0720 suggest that the magnetic 
field was not ever strong enough to decay significantly 

MM. 

We are then left with cooling. The X-ray emission 
from the INS is consistent (given the many uncertain- 
ties) with that expected for standard cooling [e.g., fnll . 
and spectral modeling finds no detectable evidence for 
non-thermal emission. This is largely consistent with the 



observed E and X-ray luminosity, as discussed above. In- 
deed, for most young neutron stars the X-ray emission 
that we see is some combination of cooling and rota- 
tion [e.g.,|72[], with the relative fractions varying depend- 
ing on source age, E, and distance (for the more distant 
sources, interstellar absorption tends to remove traces of 
soft thermal emission first). Zane et al. [69] and Kaplan 
et al. fnl suggested that RX J0720 (and the rest of the 
INS) were in fact rotation-powered pulsars either without 
radio emission or where the radio beam does not cross 
our line of sight (consistent with the very narrow beams 
found for long-period pulsars [[73Tl75ll ). but where we do 
not actually see any of the non-thermal emission. This 
is consistent with the interpretation of the Ha nebula 
around RX J 1856 as a bow-shock, which requires an en- 
ergetic particle wind to come from the neutron star (pre- 
sumably as a result of spin-down). The small deviations 
of RX J0720 from a Rayleigh-Jeans spectrum in the op- 
tical may be a hint of a non-thermal power-law peaking 
out, or they may just reflect our lack of understanding 
of the intrinsic spectrum. The spin characteristics of the 
INS are unusual compared to the bulk of the pulsar pop- 
ulation but are not entirely unheard of — we will return 
to how the INS relate to the pulsar population later. 

THE UTILITY OF THE INS 

Accepting that the INS are nearby, cooling neutron stars 
with moderate (~ 10 13 G) magnetic fields and cold (kT < 
100 eV) surfaces, the sources take on importance be- 
yond their contribution to the local neutron star popula- 
tion. We can use the INS as "physics laboratories," prob- 
ing regimes of physics not accessible experimentally (or 
sometimes even theoretically ||76|D . There are three areas 
where we can hope to derive meaningful physical con- 
straints from or explore novel regimes with the INS: (1) 
neutron star radius measurements; (2) neutron star cool- 
ing; and (3) strong magnetic fields. While measurements 
of other neutron stars can be used to address these is- 
sues, often providing complementary constraints [many 
explored in these proceedings], the INS are particularly 
well-suited. This is because they are nearby, young, and 
have no detectable non-thermal emission. Therefore they 
are bright, can still constrain cooling models, and have 
emission that is less susceptible to arbitrary decomposi- 
tion than many rotation-powered pulsars. 

The first two physical constraints come about because 
of our relative ignorance about the details of the equation 
of state (EoS) of matter at super-nuclear densities and for 
neutron-rich systems, such as those one finds in the inte- 
riors of neutron stars. We do not know some of the basic 
properties of such matter such as the constituent parti- 
cles, and options ranging from basic nuclear matter to 
pion or kaon condensates to color-superconducting quark 



matter are all possible. With each of those possibilities 
comes a range of predictions for the overall mass and ra- 
dius of the neutron star, but one finds that the radius is 
largely independent of the mass (at least for masses near 
the canonical value of 1 .4M Q ), so that a radius measure- 
ment has the hope of constraining the physical models 
of the neutron star interior IT761 p77ll . At the same time, 
the detailed predictions for the interior lead to different 
microphysical processes that affect the overall cooling of 
the neutron star, so knowledge of how neutron stars cool 
can also be used to constrain the interior IT7lll78ll . 

To constrain the radii or cooling properties of the INS, 
there are a number of observational quantities that we 
need to understand. For cooling, the goal is to place the 
INS on a plot of luminosity versus agqj. The ages can be 
estimated through two methods: either tracing the neu- 
tron stars back to likely birth locations, or through stan- 
dard pulsar spin-down assumptions. While we would ex- 
pect both of those methods to have uncertainties, the 
first is less likely to be systematically wrong. Indeed, 
the spin-down ages for RX J0720 and RX 1308 are 1- 
2 Myr, compared to kinematic ages of < 0.7 Myr mea- 
sured for those and two other sources. This could be a co- 
incidence, but it does mean that we should be cautious in 
using the spin-down ages (also see, e.g., [80]). The lumi- 
nosity needs both a distance (measured through astrom- 
etry or estimated from absorption column densities) and 
an accurate flux measurement. However, as discussed in 
Page et al. |79], the influence of the neutron star's enve- 
lope on the emergent spectrum and the cooling behav- 
ior is such that we also need to understand the elemental 
abundances, surface temperature distribution, and mag- 
netic field distribution over the surface (also see iTtsL I81L 
Yakovlev, these proceedings]). This last item has been 
the most difficult, and we will return to it below. Radius 
measurements need essentially the same data, although 
we do not need an age. But accurate distances and flux 
measurements are just as important, and again our igno- 
rance of the surface composition etc. limit any possible 
conclusions 182118311 . 

Investigating magnetic fields is a more tractable prob- 
lem. Instead of the Quantum Chromodynamics (QCD) in 
the neutron star interior, where the current theory leaves 
many areas unknown [84], the magnetic fields manifest 
Quantum Electrodynamics (QED) which is difficult but 
solvable in this regime. Here it is not so much that we 
use the neutron stars to constrain physics as we have 
two-way feedback between the theory and the observa- 
tions. The physical phase of the surface material can be 
strongly affected by the magnetic field B85I1 . and the mag- 



We could also do temperature versus age, although there are substan- 
tial difficulties in converting reliably between effective temperature and 
observed temperature [79]. 
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FIGURE 3. Energy versus magnetic field for electron and 
proton cyclotron, ground state to tightly bound (tb) or weakly 
bound (wb) states in neutral hydrogen, hydrogen ionization, 
and molecular H2 dissociation [87-89]. For fields above ~ 
10 14 G (hatched region to the right), features may partly be 
washed out due to the effects of vacuum resonance mode con- 
version 1 90]. The hatched bands show the absorption feature 
energies (corrected for gravitational redshift) for the INS (Ta- 
ble [T}. We argue 1 20] that the features in most sources are due 
to to neutral hydrogen (at least in part), the only exception be- 
ing RX J2143, which has the highest energy, and for which 
proton cyclotron absorption appears more likely. Qualitatively, 
the field strengths determined by timing for RX J0720 and 
RX J1308, indicated by the vertical lines, are consistent with 
our identifications. 



netic field also affects the radiation propagating through 
it. The vacuum is polarized, and density gradients can 
lead to transitions between the polarization states that 
have significant effects on the total spectrum [86]. Again, 
these are areas where the fundamentals are known, but 
there are no terrestrial tests to the theory in these regimes 
and hence observations of neutron stars form the best ap- 
plication of these effects. 

OPEN QUESTIONS 
Understanding the Atmospheres 

Initial attempts to model the X-ray spectra of the INS 
focused on blackbodies f56l l9lll . as they gave adequate 
fits and were simple. Including the optical/UV data it was 
clear that a single blackbody would not fit the data. In- 
stead, two blackbodies were used, where a hot, small sur- 
face fit the X-ra y da ta and a cooler, large surface fit the 
optical/UV Il2ll I92I1 . However, detailed applications of 
these models did not work, especially when one consid- 
ered the phase-dependent spectral evolution (HHHH], and 



such models are not physical. 

Instead we must consider more realistic models, tak- 
ing into account our latest knowledge of the X-ray spec- 
tra plus any additional constraints such as dipole mag- 
netic fields from X-ray timing. In interpreting the spec- 
tra, a major uncertainty is the composition. For a single 
source, this may be difficult to determine uniquely, but 
one can hope to make progress by treating the INS as 
an ensemble: ideally, it should be possible to understand 
the features (or lack thereof) in all INS with a single 
composition, appealing only to differences in tempera- 
ture and magnetic field strength (constrained by observa- 
tions where possible), which might lead to different ion- 
ization states being dominant, and possibly the formation 
of molecules or even a condensate 1871 19411 . 

There are a range of possibilities to explain the over- 
all spectra and the X-ray absorption features of the INS. 
For the absorption, the simplest models are either pro- 
ton cyclotron resonances or bound states of hydrogen, al- 
though other species are of course possible (though pos- 
sibly more difficult to explain, as gaseous atmospheres 
composed of heavier elements appear to be excluded by 
the lack of large numbers of features, and any hydrogen 
should float to the surface). 

We l20ll have attempted to determine which species 
and magnetic field is dominant for each of the INS, 
based on the energies, strengths, and widths of the ab- 
sorption features (Fig. [3]). We can find a moderately self- 
consistent picture, with both proton cyclotron and hy- 
drogen contributing. This picture agrees reasonably well 
with the magnetic fields inferred from X-ray timing al- 
though it was conceived independently. However, it is 
incomplete in at least two aspects. First, the evidence for 
multiple absorption lines Iu9l l36Tl does not fit into the 
scenario easily. Second, the evolution in the spectrum of 
RX J0720 is also not understood. See van Kerkwijk and 
Kaplan [20] for more discussion of both of these aspects, 
as well as more general issues. 

In order to fit both the X-ray and optical data simul- 
taneously, Motch et al. rf22tl and Turolla et al. l95ll con- 
sidered models that have condensed, blackbody-like sur- 
faces, where Motch et al. [22] also includes a thin layer 
of hydrogen that is optically thin at X-ray energies but 
not at optical. Getting a condensed surface is possible for 
this range of magnetic field and temperature j85ll96Tl . al- 
though it depends on the composition and the magnetic 
field distribution. Including the thin layer of hydrogen 
and varying the magnetic field allow for tuning of the op- 
tical excess to match observations. Indeed, Ho et al. 118311 
have updated this model to include partially-ionized hy- 
drogen and found a reasonable fit to the phase-averaged 
spectrum of RX J 1856 as well as the amplitude and shape 
of its pulsations [97], making this a promising avenue of 
investigation. However, we must still consider a number 
of issues addressing the physical reality of these models: 



How did such a thin layer form? Is it stable? Will it per- 
sist? Why are the layers of roughly the same thicknesses 
for different sources? 



Evolution of RX J0720 

While the initial observations of RX J0720 showed a 
smooth, blackbody-like X-ray spectrum ll2ll 19111 . later 
observations showed evidence for a phase-dependent 
absorption feature 1291 . This discrepancy was first as- 
cribed to increased sensitivity and better calibration of 
later observations, but de Vries et al. j98tl realized that 
RX J0720 was actually evolving over the course of 
months. The spectrum was changing, with the black- 
body getting harder at the same time as low-energy ab- 
sorption appeared, such that the flux was relatively con- 
stant. Simultaneously, the pulse profile evolved to in- 
crease the pulsed fraction l99ll . The change happened 
over the course of several years, but the majority of the 
variation occurred over the time span of a few months, 
from May to October of 2003. 

These observations were very puzzling, as RX J0720 
was not expected to vary on such time scales. Two sets 
of models were considered to explain the data 19811 : 
either the intrinsic spectrum of RX J0720 changed, or 
the angle at which we view the object change. They 
conclude that the second model is more likely, as no 
obvious physical motivation could be found for the first. 
Instead, free precession jlOOll could be operating. Haberl 
et al. HI in find further support for this model by fitting 
additional data and observing that the spectral changes 
seem to be reversing themselves after several years; they 
suggest a timescale of ~ 7 yr for the precession, although 
note that this was with only ~ 5 yr of spectral data, 
so the timescale was not well constrained. If true, this 
would have important implications for our study of the 
superfluid states in the object's interior [ 100]. 

However, other models are still possible. We have in- 
vestigated the spectral data and the timing data together 
i25ll . and find evidence that the changes in the spectrum 
appear more impulsive than periodic and were accompa- 
nied by a simultaneous jump in the spin frequency, with 
a fractional increase of 5 x 10~ 8 . Instead of precession, 
we interpret this as evidence for an sudden change on 
the neutron star surface accompanied by a simultaneous 
torque, followed by a slow relaxation back to the original 
spectrum. We have considered a number of mechanisms, 
from typical pulsar "glitches" that would release energy 
in the neutron star interior, to reconfiguration of the mag- 
netic field leading to changes in both the spin-down and 
in the temperature/emergent spectrum of the surface, to 
accretion of small amounts of material onto the neutron 
star surface that would both torque the star and change 



the composition. No mechanism is fully satisfactory, but 
we are attempting to include constraints from modeling 
the phase-dependence of the spectrum to try to under- 
stand the situation better [Mori et al., in prep]. We are 
also looking for similar behavior in other sources, but so 
far they seem stable [Airhart et al., these proceedings]. 



Relation to the Pulsar Population 

If we consider that the INS are moderately young 
neutron stars with ~ 10 13 G magnetic fields, we must 
ask what separates them from the bulk of the rotation- 
powered pulsar population and what their relation to that 
population is. To gain insight, we can compare with two 
groups of objects. First, the other pulsars found with 
comparable count rates (> 0.05 s _1 ) in the RASS, and 
second, the known pulsars with B > 10 13 G. 

As highlighted by its use is discovering the INS, 
the RASS is an efficient and relatively unbiased way 
to find young neutron stars. All neutron stars should 
shine in soft X-rays, either via cooling radiation or 
non-thermal processes. Taking only the thermal emis- 
sion, which should persist as a baseline, the effects 
of beaming are quite minor (compared to radio pul- 
sar searches). Aside from the 7 INS, there are 8 other 
Galact ic pu lsars with comparable count-rates [based on 
l60ll6ll[T02ll . These include many of the best-known pul- 
sars: the Crab pulsar, the Vela pulsar, PSR B0656+14, 
Geminga, and PSRs 1 055 -52, J0437 -47 1 5, B 1 95 1 +32, 
and J0538+2817 (in order of decreasing count-rate; see 
squares in Fig. |2j. Of these, the Crab is much younger 
and more energetic than the others but is significantly 
more distant, and PSR J0437— 4715 is an old recycled 
pulsar. The remaining objects are moderately young and 
nearby, much like the INS. However, the INS all have 
periods > 3 s (known for 6 of 7), while the other pul- 
sars all have periods < 0.4 s. Not only that, but the other 
pulsars all have detectable non-thermal X-ray emission. 
In fact, if we look at the characteristics that define the 
INS — young, nearby, X-ray bright, long periods — we 
cannot find any other pulsars in the ATNF catalog i65ll 
that match them. The INS are obviously then a significant 
subpopulation of objects, representing as many as half of 
the observed neutron stars within ~ 500 pc. The number 
of detected INS compared to the number of young pul- 
sars in the same volume implies a very large total number 
of similar sources in the Galaxy. There have been a num- 
ber of population estimates [e.g., 49, 66i] that try (with 
moderate success) to explore this quantitatively, but they 
are limited by the small number of sources. However, we 
must still explain why there are so many INS with simi- 
lar periods (and presumably ages and magnetic fields) in 
such a small volume. The idea of Heyl and Kulkarni 116711 . 



that magnetic-field decay could lead to a preponderance 
of nearby old magnetars, is certainly worth examining. 
While as discussed above the INS do not seem to have 
been true magnetars in the past, these conclusions used 
simple models for field decay (something still not under- 
stood in detail), and they assumed that the internal field 
was comparable to the dipolar component. If there are 
significant toroidal components to the field, these may be 
strong enough to decay but would not affect spin-down. 
Therefore, the INS could have some contribution to their 
X-ray luminosity from magnetic field decay, or magnetic 
fields could alter standard cooling behavior if they are 
strong enough to influence the heat conduction through 
to the surface ll8~Tll . 

If we the look further afield for objects that more 
closely resemble the INS, we find the so-called high- 
B radio pulsars (HBPSRs). These sources are largely 
emerging in the last seven years and their population is 
probably not complete, but we see examples with d ipola r 
fields of several xlO 13 G, almost up to 10 14 G IU03I - 
Il05ll . Some of the magnetic fields are actually higher 
than those of the INS, but the distribution of sources with 
respect to B^ip (Fig.EJ seems mostly continuous leading 
up to and past RX J0720 and RX 1308. The HBPSRs are 
younger than the INS, but tend to be much more distant 
and have higher values of E/X-my luminosity (when 
detected). This last fact could possibly be explained just 
through the usual models for pulsar evolution: assuming 
a constant magnetic field, then E evolves as B^t ~ 2 , so 

the difference of a factor of roughly 10 2 in age would 
correspond to 10 4 in E. This would then give values 
consistent with those for the INS. At the same time, as E 
dropped the dominant source of X-ray emission became 
the residual thermal emission that we see from the INS 
instead of the power-law emission that we see from the 
HBPSRs (and from other active pulsars). 

The INS may then represent a population of evolved 
HBPSRs 11691 170 ] . Previous analyses of pulsar popula- 
tions have usually not requi red p ulsars beyond magnetic 
fields of 10 13 G or so fe.g.. ll06l - H08ll . but it now seems 
apparent that the true distribution extends further in sig- 
nificant numbers (a conclusion also becoming apparent 
just from recent pulsar surveys; 01 0]). 

As discussed above, there are no detections of radio 
emission (pulsed or continuous) for the INS. The flux 
limits are reasonably low, and coupled with the small 
distances to the INS the luminosity limits are orders of 
magnitude below the luminosities of the faintest know n 
radio pulsars (such as PSR J0205+6449 in 3C 58; ifToj ). 
The non-detection of radio emission from the INS could 
be explained by the very narrow radio b eams found for 
long-period pulsars/HBPSRs (< 1% HEqI); this has 
led to large uncertainties in the predicted number of 
long -period sources [e.g., llQSh . At the same time, the 



rapid evolution of the HBPSRs across the P — P plane 
may drive them across the "death line" and terminate 
radio activity quickly, so that the objects we see may 
no longer be radio sources at all. We also look to the 
sporadically emitting RRATs, some of which seem to 
have similar spin characteristics to the INS and which 
may be a related population ll66ll . 

DISCOVERING NEW SOURCES 

With such a limited sample, the INS are invaluable for 
individual investigations but have only limited use as a 
population. In addition, each object has its own pecu- 
liarities and pathologies, which make comparison with 
the rest of the sample difficult. For instance, of the two 
brightest objects, RX J 1856 lacks an X-ray absorption 
feature but has an Ha nebula, while RX J0720 varies. 

To help remedy this situation, for several years inten- 
sive efforts have been underway to identify new isolated 
neutron stars l47l[TTo[ Pires et al. in these proceedings]. 
This is a difficult and time-consuming task, given the 
need to cross-identify X-ray sources (which often have 
poorly known positions) with optical catalogs, rej ectin g 
those sources that have plausible counterparts [e.g. Jl 1 ill . 

Currently, the majority of these efforts have yet to pro- 
duce significant results, but the first few good candidates 
are emerging. Using the SWIFT satellite, Fox ll 12ll is 
improving the positions of ROSAT sources, which facili- 
tates the search. Based on this, Rutledge et al. ll 13ll have 
identified a source which appears to be a neutron star, 
although not necessarily of the same type (young, ther- 
mally emitting) as the INS. Pires et al. [these proceed- 
ings] may have identified another from XMM-Newton 
data, although given how faint and absorbed that object 
is likely to be, confirming it as a neutron star with deep 
optical observations will be very difficult, and X-ray data 
will only yield limited information. Even so, as more of 
these objects come to light they can help us construct 
more reliable population estimates for the INS and estab- 
lish their relation to the greater neutron star population. 
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